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ABSTRACT 

Merging galaxy clusters have become one of the most important probes of dark matter, providing 
evidence for dark matter over modified gravity and even constraints on the dark matter self-interaction 
cross-section. To properly constrain the dark matter cross-section it is necessary to understand the 
dynamics of the merger, as the inferred cross-section is a function of both the velocity of the collision 
and the observed time since collision. While the best understanding of merging system dynamics 
comes from N-body simulations, these are computationally intensive and often explore only a limited 
volume of the merger phase space allowed by observed parameter uncertainty. Simple analytic models 
exist but the assumptions of these methods invalidate their results near the collision time, plus error 
propagation of the highly correlated merger parameters is unfeasible. To address these weaknesses 
I develop a Monte Carlo method to discern the merger dynamics of systems and propagate the un- 
certainty of the measured cluster parameters in an accurate and Bayesian manner. I introduce this 
method, verify it against existing numerical simulations, and apply it to two known dissociative merg- 
ers: lES 0657-558 (Bullet Cluster) and DLSCL J0916.2+2951 (Musket Ball Cluster). I find that this 
method surpasses existing analytic models, providing dynamic parameter and uncertainty estimates 
throughout the merger history, and is in better than 10% agreement with N-body simulations. This 
coupled with minimal required a priori information (subcluster mass, redshift, and projected separa- 
tion) and relatively fast computation (~6 CPU hours) makes this method ideal for large samples of 
dissociative merging clusters. 

Subject headings: galaxies: clusters: individual (lES 0657-558), galaxies: clusters: individual (DLSCL 
J0916. 2-1-2951), gravitation, methods: analytical, methods: statistical 



1. INTRODUCTION 

Merging galaxy clusters have become important 
astrophysical probes providing constraints on the 
dark matter (DM) s elf-int er action cross- s ection 
(otim: [Markevitch et all 120041: iRandah et all 120081: 
IMerten et al.l 120111: IDawson et a l."2012'). and the large- 
scale matter-antimatter ratio IStciginan 20 08i). They 
are a suspected source of e xtremely energetic cosmic 
rays (|van Weeren et al1l2010D . and the merger event po- 
tentially affects the evolution of the cluster galaxies (e.g. 
Poggianti et"all 12004 iHwang fc Lei l2009t IChung et al.l 
2009f ). All of the respective astrophysical conclusions 
drawn from merging clusters depend on the specific 
dynamic properties of a given merger. 

For example, the subclass of dissociative mergers, in 
which the collisional cluster gas has become dissoci- 
ated from the near collisionless galaxies and dark mat- 
ter, provides four ways of constrai ning the dark mat- 
ter self-interaction c ross-section (Markevitch et al . 2004: 
IRandall et 311120081 ). The best constraints come from 
studying the mass-to-light ratios {M/L) of the subclus- 
terfl, and the ofi :set between the collisionless galaxies 
and dark matter ([Markevitch et al.l 120041 : IRandall et al.l 
|2008[ ). Both constraints directly depend on the merger 
dynamics. First the relative collision velocity will af- 
fect the expected momentum transfer between each sub- 

|wadawson@ucdavis.edu 

^ I define subcluster as either one of the two colliding clusters, ir- 
respective of mass, and I define cluster as the whole two-subcluster 
system. 



cluster's dark matter particles which will in turn af- 
fect the expected dark matter mass transfer from the 
smaller subcluster to the larger subcluster ultimately 
affec ting the expected mass to light ratios of the clus- 
ters (|Markevitch et al.ll2004[ ). Second the observed time- 
since-collisior0 will affect the expected galaxy-dark mat- 
ter offset, with the offset increasing with time (for ctdm > 
0) as the system approaches the apoapsis (or turn-around 
distance), dmax- Additionally it is important to know 
the velocity so that dark matter ca ndidates wi t h ve- 
locity dependent cross- sections (e.g. IColin et "aTl 120021 : 
IVogelsberger et al]l2012[ ) can be constrained. 

However there is no way to directly observe the dy- 
namic merger parameters of principal interest: the three- 
dimensional relative velocity (wsd) and separation (dan) 
of the subclusters as a function of time, their maximum 
separation (dmax), the period between collisions (T), and 
the time-since-collision [TSC). Observations are gener- 
ally limited to: their projected separation (rfproj), the 
line-of-sight (LOS) velocity of each subcluster {vt) as in- 
ferred from their redshifts, and their mass {Mi) or pro- 
jected surface mass density profile. In addition to the ob- 
vious inability to measure a change in the merger state, it 
is difficult to constrain the dynamic parameters of inter- 
est even in the observed state. This is due to the general 
inability to constrain the angle of the merger axis with 
respect to the plane of the sky (a), see Figure [1] 

For the Bullet Cluster it was originally thought that 

^ I define the time of collision to be the time of the first pericen- 
tric passage. 
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estimates of the Mach number of the cluster merger 
throuRh X-ray ob servations of the gas shock feature (e.g. 
iMarkevitchI 120061 ) could be used to estimate wsd and, 
in conjunction with measurements o f the relative LOS 
velocities, then estimate a. However, iSpringel fc Farra^ 
12007) showed that the Mach number only translates to 
an upper limit on i;3d, and in the case of the Bullet Clus- 
ter they showed that the Mach inferred velocity could be 
a factor of ^ 2 larger than the true wsd . There is poten- 
tial for constraining a using polarization measurements 
of radio relics ( Ensslin e t al. 1998], which arc associated 
with some cluste r mergers (e.g . Ivan Weeren et al.ll20lll 
but not all (e.g. iRussell et al.|[2^11f) . Even if for some 
mergers radio relics provide constraints on a, dynamic 
models are still needed in order to ascertain the dynamic 
properties of the merger throughout time. 

The two most prevalent methods for ascertaining the 
dynamics of observed merging systems are the timing 
argument and N-body simulations. The timing argu- 
ment is based on the solution to the equations of mo- 
tion of two gravitating point masses, with the cosmo- 
logical constraint that as 2; — > oo the separation of the 
two masses dgn (for an exposition of this method 
see jPecblcs 19 931) . The tim ing argument was first used 
by iKahn fc Woltieii (|19590 to study the system of the 
Milky Way and M31, and firs t applied to binary cluster 
systems bv lBeers et al.l ()1982l ). It has recently been ap- 
plied to s everal dissociative m ergers, includ ing the Bullet 
Cluster (Barrcna ct al. '200 2D, Abell 5 20 (jCirardi et al.l 
[20081 Abell 21 63 ( gourdin et al.ll20lil ). and Abeh 1758N 
(jBoschin et al.ll2012D . N-body simulations of observed 
dissociative mergers have been limited to the Bullet Clus- 
ter (lES 065 7-558) and have come in two variants: hy- 
drodynamic ("Springel fc Farrar"2007:'Milosav lievic et ahl 
[200 7: Mastropictro & Burkcrt 2008), and self interacting 
dark matter (S IDM) plus collisionless galaxy particles 
([Randall et al.ll2008. ). 

While the timing argument method is easy to use, its 
inherent assumptions result in non-negligible error for 
dissociative systems. Most importantly the timing ar- 
gument method assumes two point mass particles; this 
assumption begins to break down as the two subclusters 
overlap and results in divergent solutions as the subclus- 
ters near collision. Since most dissociative mergers are 
observed with the two-subcluster halos overlapping and 
dark matter constraints depend on the merger dynam- 
ics near collision, application of the traditional timing 
argument to dissociative merger s is lim ited and should 
be done with caution. iNusseii (|2008l ) addressed this 
weakness of the traditional timing argume n t by s ubsti- 
tuting truncated iNavarro. Frenk. fc Whit d (| 19961 . here- 
after NEW) halos and numerically solving the equa- 
tions of motion. Another weakness of the timing argu- 
ment is that its main constraint requires the assumption 
that the subcluster masses a re constant since the begin- 
ning of the universe. While lAneus fc McGaughl ([20081 ) 
have noted this problem with the initial conditions of 
N-body simulations and proposed a solution based on 
esti mating the mass accr etion histories of the clusters 
(e.g. IWechsler et "al1l2002| ) . their correction is incompati- 
ble with the timing argument method as this would add 
a second differential term to the equations of motion. Fi- 
nally, the large covariance between the merger parame- 
ters plus the complexity of the equations of motion makes 




Fig. 1. — The generic two-halo merger configuration assumed in 
this work. Observable parameters are shown in dark blue, and in- 
clude the mass of each halo (Mi), the projected separation (rfproj), 
and the line of sight (LOS) velocity components {vi) as deter- 
mined from the halo redshifts. The generally unknown parame- 
ters of the mergers are shown in light blue, and include the angle 
of the merger axis with respect to the plane of the sky (a), and 
the three-dimensional separation (oiaD) and velocity components 
(1)30. ). Note that while just the outgoing scenario is shown in this 
figure, the method also considers the incoming scenario. 

propagation of errors in the timing argument formalism 
untenable. This has resulted in a lack of certainty with 
timing argument results, leaving most users to run a few 
scenarios in an effo rt to roug hly bound the range of pos- 
sible solutions (e.g. iBoschin et al. 2012|). 

N-body simulations provide the most accurate descrip- 
tion of merger dynamics, however they are computation- 
ally expensive which results in their application being 
limited. Despite eleven currently confirmed dissociative 
merger^! only the Bullet Cluster has been modeled with 
N-body simulations, whereas most of these have been 
analyzed with the timing argument method. Existing 
N-body simulation strategies to ascertain the dynamic 
properties of mergers are incapable of keeping up with 
the current faster than exponential rate of discovery. 
Even for the case of the Bullet Cluster the N-body anal- 
yses have been limited as far as mapping out the merger 
dynamic phase space allowed by the uncertainty of the 
obse rvations, with at most 15 diff e rent scenarios being 
run (jMastropietro fc Burkeri2008D . iGomez et all (|2012[ ) 
have come the closest to addressing this issue in their in- 
vestigation of potential dissociative mergers (A665 and 
AS1063) through the use of simpl ified scale-free num eri- 
cal simulations of the mergers (see IGomez eTaII[2000l for 
details). However, they have still had to severely limit 
the phase space probed (fixing merger parameters such 
as the initial relative velocity and subcluster-subclustcr 
mass ratio); thus admittedly this approach enables con- 

3 (1) Bull et Cluster ijClowe et al.l [2003 ): (2) A520 

JMahdavi et al.| [20071); (3) MAC S J0025.4-1222 l IBradac et all 
"2008); (4) A1240 (Barrena et al. 200l) : (5) ZwC L 0008.8-1-5215 
(van Weeren et al. 2011); (6) A2744 JMerten et al. 201 lj); (7 ) 
A2163 t Okabe et al. , 2011); (8) A1758N (Rago zzine et aTTISoTa) ; 
(9) Musket Ball Clu ster (D awson et al. 2012); (10) ACT-CL 
J01024915 ( Menantea u et al.l I2012D : (11) MACS J0717.5-I-3745 
HMroczkowski et al. 2 01^ ) 
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struction of plausible models, but not a thorough ac- 
counting of possible or likely models. 

With these weaknesses in mind I present a new 
methocfl for analyzing the dynamics of observed disso- 
ciative mergers. My primary objectives are to 1) obtain 
a solution valid near the collision state, 2) fully estimate 
the covariance matrix for the merger parameters, 3) be 
able to analyze a dissociative merger on the order of a day 
using a typical desktop computer, and 4) obtain better 
than 10% accuracy; all assuming that only the most gen- 
eral merger observables and their uncertainty are known: 
mass of each subcluster, redshift of each subcluster, and 
projected separation of the subclusters. 

In ij2] I define a method for analyzing the dynamics 
of observed dissociative mergers. In i )2.3l I verify this 
method with existing results from hydrodynamic N-body 
simulations. In ij3] I apply this method to the Bullet 
Cluster and in f|4] I apply this method to the Musket 
Ball Cluster (DLSCL J0816.2-K2951) and contrast its 
dynamics with those of the Bullet Cluster. Finally in 
Sj5] I summarize my findings, discuss their implications 
for the constraints on dark matter and suggest other 
science that will benefit from the introduced method. 
Throughout this paper I assume JIa = 0.7, flm = 0.3, 
and H = 70kms~i Mpc"^ 

2. METHOD 

In order to obtain a valid solution of the system dynam- 
ics near the collision state I use a model of two spherically 
symmetric NFW halos, rather than point masses. I incor- 
porate this model in a standard Monte Carlo implemen- 
tation: draw randomly from the observables' probability 
density functions (PDF's) to generate a possible realiza- 
tion of the merger, use the model to calculate merger 
properties of interest, apply multiple priors, store these 
likelihood weighted results as a representative random 
draw of their PDF, and repeat. The final result is a mul- 
tidimensional PDF for the dynamic parameters of the 
merger. This method agrees well with hydrodynamic 
simulations, tJ2.3[ and satisfies the speed and accuracy 
objectives outlined in the Introduction. 

2.1. Model 

The general basis of the model is a collisionless two 
body system with the mass of each body mutually con- 
served throughout the merger. The model requires min- 
imal input: the mass of each subcluster, the redshift of 
each subcluster, and the projected separation of the sub- 
clusters (along with associated uncertainties). It assumes 
conservation of energy and zero angular momentum. The 
model also assumes that the maximum relative velocity 
of the two bodies is the free-fall velocity of the system 
assuming their observed mass. In the remainder of this 
subsection I will discuss in detail these general assump- 
tions, their justification, and their implications. 

I model the system using two spherically symmetric 
NFW halos truncated at r2ocEl. By default the concen- 
tration of each halo is determined by the halo's mass 

This method is similar to the one used bv lDawson et ahl II2012I ). 
although with several improvements (see iffl . 

^ r200 is defined as the radius of the spherical region within 
which the average density is 200 times the critical density at the 
respective redshift. 



via th e mass-concentration scaling relation of lDufiv et al.l 
(|2008f) . This is not a requirement of the model though, 
and measured concentrations can be used, as in the case 
of iJ3.1l The dynamic parameter results are relatively 
insensitive to the assumed concentration of the subclus- 
ters. Take for example the case of ij2.3l with user speci- 
fied concentration s of Ci = 1.94 and C2 — 7.12: if instead 
iDuffv et"al] ()2008D inferred concentrations Ci = 3.44 and 
C2 = 2.75 (~ 200% difference for both) are used, the dif- 
ference in the estimated fsoCicoi) and TSC are both less 
than 6%. 

The model assumes that the mass of each subcluster 
is constant and equal to the observed masfl While this 
assumption is also used in the timing argument method, 
it is more reasonable for this method since the bulk of the 
results are calculated between the observed state and the 
collision state, typically lasting < 1 Gyr. This is about an 
order of magnitude shorter than the typical timescales of 
the timing argument method thus the new method is less 
susceptible to error due to neglecting growth of structure. 

The model assumes that the energy of the two-halo 
system is conserved, and consists only of their mutual 
kinetic and potential energies. The kinetic energy of the 
system is K{t) — O.bfiv^nit)^ , where /i is the reduced 
mass of the system and vsnit) is the relative physical ve- 
locity of the two subclusters at time t. The potential en- 
ergy of the system is assumed to be purely gravitational 
and is derived in Appendix [Xj Since the model assumes 
zero impact parameter there is no ro tational kinetic en- 
ergy term. iMastropietro fc Burkertl (.2008.) find that a 
moderate impact parameter of ~ 0.1r2oo has less than a 
1% effect on the merger velocity, thus this assumption 
should have negligible effect for the case of dissociative 
mergers which must have had relatively small impact pa- 
rameters in order to dissociate the bulk of their gas. 

For the relative velocity of the two subclusters I apply a 
fiat prior from zero to the free-fall velocity of the subclus- 
ters, assuming their observed mass. This will result in an 
overestimate of the maximum possible relative velocity, 
due to the neglect of mass accretion. It is conceivable 
that this prior could be tightened using the maximum 
relative velocities observed in cosmological N-body sim- 
ulations as a function of subcluster masses and redshift. 
Another possibility for tightening the prior would be to 
analytically estimate the free-fall velocity acco unting for 
mass accretion (e.g. lAngus &: McGaughl I2008D . An ad- 
vantage of the Monte Carlo approach taken with this 
method is that additional priors can be applied as more 
knowledge becomes available without the need to rerun 
the analysis, so I opt for a conservative default approach. 

The model ignores the effects of surrounding large scale 
structu re and simply treats the two-body system. As 
'Nusserl (|2008t ) shows, a global overdense region (10 times 
denser than the background) engulfing the system only 
affects the dynamics substantially for extreme collision 
velocities (~ 4500 km s""'^). While global overdensities 
may be disregarded it is not clear that the effects of 
nearby structures can be disregarded, e.g. as in the case 
three body systems. Thus this method should be applied 
with caution to complex cluster mergers. 

The model also ignores dynamical friction. 

^ For subcluster mass I refer to M20O1 which is the mass of the 
individual subcluster enclosed within a radius of r200 ■ 
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iFarrar fc Rosed ()2007D found that including dynamical 
friction accounted for an ~10% reduction in the inferred 
collision velocity of the Bullet Cluster in their analytic 
treatment. This is potentially concerning since dynam- 
ical friction is inversely proportional to the relative 
velocity of the merger, thus it may become even more 
important for mergers slower than the Bullet Cluser. 
However in ^2.3\ I compare the results of my method 
with those from hydrodynamic N-body simulations and 
show that the net effect of all simplifications (including 
ignoring dynamical friction, tidal stripping of mass 
and gas mass lost during the collision) are negligible, 
suggesting that dynami c friction is less i mport ant than 
the analytic estimates of lFarrar fc RoseiJ ()2007[ ) suggest. 

2.2. Monte Carlo Analysis 

In this section I discuss the details of the Monte Carlo 
analysis workflow. I chose to implement a Monte Carlo 
analysis because the high degree of correlation among 
the many merger dynamic parameters made traditional 
propagation of errors unfeasible. A Monte Carlo analysis 
has the added advantage of easily enabling application 
of different combinations of priors ex post facto, see e.g. 

The analysis begins by randomly drawing from the 
PDF's of the merger observables: mass of each subcluster 
(-^200^7 redshift of each subcluster (zi), and projected 
separation of the subclusters (dproj)- The potential en- 
ergy, V (see Appendix \^ , at the time of the collision is 
used to calculate the maximum relative velocity, 



locity is calculated. 



V{r = 0). 



The velocity of each subcluster relative to us is esti- 
mated from its redshift. 



(1 



1 



C, 



where c is the speed of light. The relative radial velocity 
of the subclusters is calculated from their redshifts. 



Wrad(^obs) 



1^2 - Vl\ 



Since the angle of the merger axis with respect to the 
plane of the sky, a, is unconstrained without prior knowl- 
edge of the three-dimensional relative velocity, I assume 
that all merger directions are equally probable. How- 
ever, projection effects result in PDF (a) — cos(a). Due 
to symmetry it is only necessary to analyze the range 
< a < 90 degrees. I draw randomly from this PDF 
for each realization. This enables the calculation of the 
three-dimensional relative velocity in the observed state. 



W3D(iobs) = Wrad(tobs)/ sin(a). 



(1) 



as well as the observed three-dimensional separation of 
the subclusters, 



dsDitohs) = c?proj/ cos(a). 



(2) 



If W3D(^obs) > ^^3D„,ax: then this realization of the 
merger is discarded; otherwise the relative collision ve- 



^30(^001) = y W3D(iobs)^ + - [V{tohs) - V^(icol)]- 

Similarly if U3d(^coi) > ^3D,„axJ then this realization is 
discarded. 

The change in time. At, between two separations is 
given by 

n dr 

At = 



I define the time-since-coUision (TSC) as the time it 
takes the subclusters to traverse from zero separation 
to their physical separation in the observed state, dsu- 
Because there is a potential degeneracy in whether the 
subclusters are "outgoing" (approaching the apoapsis af- 
ter collision) or "incoming" (on a return trajectory after 
colliding and reaching the apoapsis); I solve for both of 
these cases, TSCq and TSCi respectively. In determin- 
ing TSCi it is useful to define the period, T, of the sys- 
tem. I define T to be the time between collisions. 



T 



dr 



^liE-V{r)) 



where dmax is the distance from zero separation to the 
apoapsis, when E — V . Thus, 

TSCi =T- TSCq. 

During the Monte Carlo analysis any realizations with 
TSCq greater than the age of the Universe at the clus- 
ter redshift are discarded. A similar flat prior is applied 
when calculating the statistics of TSCi. To this regard 
some insight into the likelihood of the system being in an 
"outgoing" or "incoming" state can be gained by calcu- 
lating the fraction of realizations with TSCi less than the 
age of the Universe at the cluster redshift. Conceivably 
these temporal priors could be strengthened, requiring 
that the time to first collision (T) plus the respective 
TSC be less than the age of the Universe at the clus- 
ter redshift, in a fashion similar to the timing argument. 
However, as with the timing argument model, the model 
of W2.1\ becomes less valid over time-scales approaching 
the age of the Universe. Thus I use the more conservative 
prior by default. 

Since the majority of the merger time is spent at large 
separations, due to lower relative velocities, observations 
of the system are more likely near apoapsis than near 
the collision. Thus the probability of each realization is 
convolved with the prior 



PDF(rS'Co) = 2 



TSCq 

T 



(3) 



There are likely selection effects which complicate this 
PDF, since it can be imagined that the X-ray luminos- 
ity is greatest near the time of the collision (see e.g. 
I Randall et al .11200^ . However this information is rarely 
if ever known, thus it is not included by default. In > j3.2.2l 
I show how additional temporal priors, based on similar 
effects, may be effectively applied to the results of the 
analysis ex post facto. 
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TABLE 1 

Comparison of the model results with 
hydrodynamic simulation results 



Work 


f3D(tcol) 


TSCo 




kms~^ 


Gyr 


Soringel & Farrar C2007~) 


3400 


0.18 


This work 


3380 


0.198 


% Error 


0.70 


10 



The end result of this method is a 13 dimensional pos- 
terior PDF of an array of cluster merger parameters, see 
for example Appendix [BJ Finally to compact the results 
I use the biweight-statistic (generally more robust and 
less sensitive to abnormally tailed distributions than the 
median or mean) and bias-corrected percent confidence 
limits (Beers et al. 1990) applied to the marginalized pa- 
rameter distributions of the valid realizations, see for ex- 
ample Table [3l 

2.3. Comparison with Hydrodynamic Simulations 

For the purposes of checking the physical assumptions 
of the model I reanalyze the iSpringel fc Farrail ()2007f ) 
model of the Bullet Cluster, comparing my dynamic pa- 
rameter estimates with their hydrodynamic N-body sim- 
ulation based estimates. For this analysis I run just their 
single case through the model (i.e. I do not perform a 
Monte Carlo analysis). Based on private communication 
from D. Clowe they represent the "main" and "bullet" 
subclusters as NFW halos with M2001 = 1-5 x IO^'^Mq, 
ci = 1.94, M2002 = 1-5 X 10" Mq, and C2 = 7.12, re- 
spectively. They note that the gas properties of their 
simulation most closely match the observed Bullet Clus- 
ter gas properties for the time step corresponding to a 
subcluster separation of o?3d = 625 kpc and relative ve- 
locity of W3D(^obs) = 2700 km s~^. I define this as the 
"observed" state and use these values as input to the 
model. I compare t he results of the two mo dels in Ta- 
ble m assuming the ISpringel fc Farrail ()2007D results as 
truth when calculating the percent error. There is bet- 
ter than 1% agreement between the collision velocities 
and 10% agreement between the TSCq. I expect better 
agreement between the collision velocity estimates, than 
between the TSCq estimates, since the velocity calcu- 
lation comes from simply comparing the observed and 
collision states of the merger whereas the TSC calcula- 
tion requires the integration of the equations of motion 
between these two states. These results essentially vali- 
date many of the simplifying assumptions of the model 
(conservation of energy, and ignoring the affects of dy- 
namical friction, tidal stripping of dark matter and gas 
during the collision). 

3. BULLET CLUSTER DYNAMICS 

The Bullet Cluster is the prime candidate for first ap- 
plication of the method as it is one of the best studied 
dissociative mergers. It has a wealth of observational 
data necessary for input to the model, as discussed in 
§3.11 plus supplementary data which enables additional 
posterior priors, as discussed in §3.2.21 

3.1. Bullet Cluster Observed System Properties 



TABLE 2 

Bullet Cluster parameter input 



Parameter 


Units 


M 


a 


Rcf. 


M2001 


10" Mq 


15 


1.5^ 


1 


ci 




7.2 


b 


1 


M2002 


10" Mq 


1.5 


0.15=^ 


1 


C2 




2.0 


b 


1 


Zl 




0.356211 


0.000342 


2 


Z2 




0.360174 


0.000204 


2 


'^proj 


kpc 


720 


25 
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References. — (1) Springel & Farrar 2007; (2) 
Barrena et al. 2002; (3) Bradac et al. 2006. 



Note. — A Gaussian distribution with mean, fi, 
and standard deviation, a, is assumed for all param- 
eters with quoted respective values. The mass, M200, 
and concentration, c, are the defining properties of as- 
sumed spherically symmetric NFW halos. 

^ Estimated to be 10%, base d one the error magn itude 
of M(< 250 kpc) reported in lBradac eFal] (12009) . 
^ No errors were presented in the reference. A sin- 
gle concentration value was used for all Monte Carlo 
realizations. 

I summarize the observed Bullet Cluster parameters 
used as input to my analysis in Table [5] The full PDF's 
of these input parameters have not been published so 
I simply assume Gaussian distributions. I refer to the 
main subcluster as halo 1 and the "bullet" subcluster 
as halo 2. For the mass and concentration of each sub- 
cluster I use the most re cently reported estimates from 
ISpringel fc FarraJ (I2007D. based upon strong and weak 
lensing estimates (jBradac et al.l i2006). However, they 
do not present errors for these quantities so for the mass 
I estimate the l-tr errors to be 10% of the mass, since this 
is approximately th e magnitude of the error reported by 
iBradac et all ([200I for M(< 250 kpc). There is no pub- 
lished estimate for the uncertainty of the concentrations 
of the NFW model fits, q, so I simply assume the con- 
centrations to be known quantities (as noted in H2.1\ the 
results are relatively insensitive to the assumed concen- 
trations). The redshifts of the main and "bullet" subclus- 
ters are esti mated from 71 and 7 s pectroscopic members, 
respectively (jBarrena et al.l 12002^. The projected sepa- 
ration of the mass peaks is determined from str ong and 
weak lensing measurements (jBradac et al.l 120061 ). and is 
essentially the same as the separation of the subclusters' 
galaxy centroids. For each Monte Carlo realization in- 
dividual values are drawn randomly from each of these 
assumed Gaussian distributions. 

3.2. Bullet Cluster System Dynamics Results 

I first analyze the Bullet Cluster with the Monte Carlo 
analysis method and default priors discussed in 321 high- 
lighting the complexity of merger dynamics and the inap- 
propriateness of analyzing a small sample of select merger 
scenarios. In i J3.2.2l I incorporate additional constraints 
provided by the observed strong X-ray shock front plus 
boosted temperature and luminosity. I discuss this prior 
information and apply it ex post facto to the default prior 
resuhs of g32H 

I perform the analysis with 2,000,000 Monte Carlo re- 
alizations. Parameter estimates converge to better than 
a fraction of a percent with only 20,000 realizations 
(~ 6 CPU hours). I run a factor of a hundred more since 
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Fig. 2. — The posterior of the Bullet Cluster's time-since- 
coUision TSCq and 1^30(^001) parameters is shown in grayscale with 
dark and light blue contours representing 68% and 95% confidence, 
respectively. The green-scale triangles are from a subsample of the 
Monte Carlo population, which jointly satisfies the requirement 
of being drawn from itO.Olcr of the mean of each of the input 
parameters, i.e. the "most likely" values for the input parame- 
ters. Despite representing the most probable input parameter val- 
ues there is considerable spread in the inferred output parameters, 
with the subsample clearly tracing the ridge of the distribution. 
The saturation of the triangles increases with increasing a, from 
15— 85 degrees. The purple-scale circles arc from a subsample near 
the bi-weight location of o = 52 ± 0.0002 degrees, with the satu- 
ration of the circles increasing with increasing M20O1 • While the 
length of the distribution is predominantly caused by uncertainty 
in a the width is predominantly caused by uncertainty in the in- 
put parameters. Despite the Bullet Cluster being one of the best 
measured dissociative mergers there is still considerable and com- 
plex uncertainty in its merger parameters, predominantly due to 
uncertainty in a. 

it was computationally inexpensive and it provides a data 
sample to which I can apply any number of conceivable 
posterior PDF's and still maintain sub-percent statistical 
accuracy. 

3.2.1. Default Priors 

The main results of this analysis are that: 1) there is a 
great degree of covariance between the geometry, veloc- 
ity, and time parameters of the merger, and 2) models 
of the system which disregard the uncertainty of a will 
catastrophically fail to capture the true uncertainty in 
the dynamic parameters. 

The two-dimensional PDF of Figure [5] exemplifies the 
complexity of the covariance between the various merger 
parametersQ. The shape of the PDF is most easily under- 
stood in terms of the parameters' dependence on a. This 
dependence is illustrated by the green-scale triangles that 
represent a subsample of the Monte Carlo population, 
which jointly satisfies the requirement of being drawn 
from ibO.Olcr of the mean of each of the input parameters, 
i.e. the "most likely" values for the input parameters. 
The saturation of the triangles increases with increasing 
a, from 15-85 degrees, clearly showing a monotonically 
increasing relationship with TSCq (see also Figure [Tl)) . 
For small a (light green triangles), Equation[T]states that 
f3D(^obs) must be large thus v^Ditcoi) must also be large. 

Similar degrees of complex covariance exist for the other ge- 
ometry, velocity and time parameters, see e.g. the results array in 
Appendix iBl 



and since Equation[5]states that daDC^obs) approaches the 
minimum possible observed separation, lipi-oj , the TSCo 
must approach a minimum. Conversely for large a (dark 
green triangles), (i3D(iobs) becomes large increasing the 
time required to reach the observed state, and despite 
W3D(^obs) approaching the minimum 'yrad(iobs) the colli- 
sion velocity must increase for the subclusters to have 
been able to reach the larger dsuitohs)- 

The bulk of the uncertainty in the geometry, velocity 
and time parameters is due to the uncertainty of a. This 
is exemplified by the fact that the green-scale triangles in 
Figure [5] closely trace the extent of the ridge line of the 
two-dimensional distribution (i.e. span the bulk of the 
uncertainty) . Conversely the "width" of the distribution 
is predominantly due to uncertainty in the input param- 
eters. This is exemplified by the purple circles of Figure 
[21 which are for a near constant a yet randomly sam- 
ple the M2001 distribution. The saturation of the circles 
increase with increasing mass. 

The inability to directly measure a, coupled with its 
strong degree of correlation with the other dynamic 
parameters, makes it the dominant source of uncer- 
tainty. While it was originally believed that the three- 
dimensional merger velocity as inferred from the X-ray 
shock feature could be coupled with the redshift deter- 
mined radial velocity to measure a, iSpringel fc Farrail 
12007') showed that the X -ray shock inferred velocity 
significantly overestimates the true three-dimensional 
merger velocity. So at best this information can weakly 
constrain a, and in the case of the Bullet Cluster the X- 
ray shock inferred velocity is significantly greater than 
the free-fall velocity, wsD^axi t hus it provides no addi- 
tional constraining power. In I discuss how the re- 
sults of this method can be used in conjunction with 
N-body simulations to limit the computational impact 
of accounting for the uncertainty in a. 

3.2.2. Added Temporal Prior 

One of the advantages of this Monte Carlo method 
is that additional constraints are easily incorporated ex 
post facto. An example of such constraints in the case of 
the Bullet Cluster is the observed X-ray shock front and 
factor of 2.4 gr eater X-ra y estimated mass to lensing es- 
timated mass (jMarkevit ch 2006), due to merger related 
X-ray temperature and luminosity boo st. Hydrodynamic 
simulations of merging clusters (e.g. iRicker fc SarazinI 
I2OOII : iRandall et al.ll2002l ) suggest that such transient ef- 
fects last of order the X-ray sound crossing time. Since 
simulations show negligible difference between the time 
scales of the two I chose to con struct a prior based o n the 
observed temperature boost. IRandall et al.l (|2002[ ) find 
that the full-width-half-max (FWHM) duration of the 
temperature boost is ~ OAtsc with the entire boost du- 
ration being ~ lAtgc, where tsc is the sound crossing time 
of the more massive of the two subclusters. The peak of 
this boost roughly coincides with the time of the colli- 
sion, as defined in 21 Given the M2001 = 15 x lO^'' M© 
and temperature Tx = 14keV of the "main" subcluster 
(|Markevitchl2006( ) . the tsc — 1 Gyr. I construct a sigmoid 
function for the TSC prior PDF based on the observed 
temperature boost. 



Dynamics of Merging Clusters 



7 



TABLE 3 

Bullet Cluster parameter estimates 



Parameter 


Units 




Default Priors 




Default + Added Temporal Priors 






Location^ 


68% LCL-UCL'> 


95% LCL-UCL'' 


Location'' 


68% LCL-UCL'' 


95% LCL-UCL'' 


M20O1 


lO'-* M0 


15.1 


13.6 - 16.6 


12.1 - 18.1 


10. z 


lO-O — iU.O 




iVl2002 




1 5 


13 — 16 


19 IS 


1.5 


1.4 - 1.7 


1.2 - 1.8 


Zl 


0.3562 


0.3559 - 0.3566 


0.3556 - 0.3569 


0.3562 


0.3559 - 0.3566 


0.3556 - 0.3569 


Z2 




0.3602 


0.36 - 0.3604 


0.3598 - 0.3606 


0.3602 


0.36 - 0.3604 


0.3598 - 0.3606 




Mpc 


0.72 


0.69 - 0.76 


0.65 - 0.80 


0.72 


0.68 - 0.75 


0.64 - 0.79 


a 


degree 


52 


31 - 72 


19 - 83 


30 


21 - 43 


16 - 56 




Mpc 


1.2 


0.9 - 2.7 


0.7 - 6.7 


0.8 


0.8 - 1.0 


0.7 - 1.3 




Mpc 


1.5 


1.2 - 2.4 


1.1 - 5.6 


1.4 


1.2 - 1.9 


1.1 - 3.5 


l'3I) (*obs ) 


km 


1100 


890 - 1700 


780 - 2600 


1800 


1400 - 2600 


1100 - 3300 




km 


3200 


3000 - 3900 


2800 - 4300 


3100 


2900 - 3500 


2800 - 3900 


TSCo 


Gyr 


0.5 


0.3 - 1.2 


0.3 - 3.8 


0.3 


0.2 - 0.4 


0.2 - 0.6 


TSCi" 


Gyr 


1.5 


1.2 - 2.4 


1.0 - 7.5 


1.4 


1.2 - 2.3 


1.0 - 5.3 


T 


Gyr 


2.0 


1.6 - 3.2 


1.4 - 11.0 


1.9 


1.6 - 2.3 


1.4 - 5.0 



^ Biweight-statistic location (see e.g.lBe ers et al.lll99g) . 

Bias-corrected lower and upper confidence limits, LCL and UCL respectively (see e.g. lBeers et al.lll99CII ). 

For the case of the Default + Added Temporal Prior, none of t he real izations have a valid TSCi, meaning that the Bullet 
Cluster is being observed in the "outgoing" state, as discussed in ^3.2.21 




0.1 
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Fig. 3. — The posterior of the Bullet Cluster's TSCo and 
D3d(4coi) parameters after application of an additional temporal 
prior based on X-ray observations of the Bullet Cluster (grayscale). 
Dark and light blue contours representing 68% and 95% confidence, 
respectively. The added temporal prior significantly improves the 
constraint on the merger parameters (compare v^ith Figure [2]|. 
The black diamond represents the Springel & Farrar (2007) hy- 
drodynamic simulation result for their defined "observed state", 
whose X-ray properties best match the observed X-ray properties; 
'i3D=625kpc for this state. The green bar shows their result for 
d-so between 700 to 900 kpc, which is more in line with the observed 
dproj = 720 ± 25 kpc (assuming < o < 35 degrees). 

where a is the FWHM of the duration of the temperature 
boost and b is the entire boost duration. I chose a sigmoid 
function over a simple st ep function sin ce the tempera- 
ture boost predicted by I Randall et al.l p002i) does not 
end abruptly. This prior is coupled with the previously 
discussed TSC prior (Equation [3]). 

Application of this prior significantly improves the un- 
certainty in TSCo (180% to 67%) and U3D(icoi) (28% 
to 19%), compare Figure [5] with Figure |31 It essen- 
tially removes the possibility of a TSC > 0.7 Gyr. As 
expected from the a dependence shown by the green tri- 
angles in Figure [5] this prior also reduces the likelihood 
of a > 50 degrees, which in turn affects both the loca- 
tion and uncertainty of d^B and f3D(iobs), see Table [3l 



The remaining parameter estimates are predominantly 
unaffected by the prior, with only a few having their 
confidence limits affected as the result of their high end 
low probability tails being down weighted. Additionally 
there is now essentially zero probability that the bullet 
subcluster has reached the apoapsis and is on a return 
trajectory, since the 95% lower confidence limit of TSCi 
is 1.0 Gyr (see Table [31) and the prior essentially removes 
the possibility of a TSC > 0.7 Gyr. I present the full re- 
sults array, which includes the default analysis prior com- 
bined with this temperature boost prior, in Appendix iBl 
and the compact parameter est imates in Table [3l 

According to th is an alysis the lSpringel fc Farrad (|2007l ) 
results stated in ij2.3l seem unlikely (see the black dia- 
mond of Figure [3]) , however this is simply due to their 
definition of the "observed state" . They define the ob- 
served state to be when their simulated X-ray properties 
most closely match the observed X-ray properties, yet the 
separation of the halos in this state is only d^u = 625 kpc; 
this is less than t he observed dproj = 720 ± 25 kpc 
(IBradac et al.l 120061) . If we instead consider their esti- 
mate of TSCo for dsD between 700 to 900 kpc (corre- 
sponding to dproj — 700 and < a < 35 degrees), then 
0.24 < TSCo < 0.33 Gyr (see green bar of Figure Ej). 
This brings their result in line with the results of this 
method, as expected by the agreement prese nted in ii2.3l 
Note that the general conclusion of Spring el fc Farran 
(|2007f ). that the shock speed greatly overestimates the 
actual relative speed of the subclusters, remains valid 
regardless of which "observed state" is used. 

4. MUSKET BALL CLUSTER DYNAMICS 

I also apply the method to the Musket Ball Cluster, 
with the objective of updating an existing analysis and 
comparing this system with the Bullet Cluster. A pre- 
liminary analy sis of the system dynamics using a similar 
method (Daw son et al.l 12012) suggested that the Mus- 
ket Ball Cluster merger is ~3-5 times further progressed 
than other confirmed dissociative mergers. However, 
that analysis treated the two merging subclusters as uni- 
form density spheres and also failed to account for the 
temporal phase-space PDF (Equation [S]) . Additionally 
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Fi g. 4. — Weak lensing mass PDF's of the Musket Ball subclus- 
ters l|Dawson et al.|[2"0ia ). 

the claim that the Musket Ball Cluster is both slower and 
further progressed than the Bullet Cluster was based on 
comparing the Musket B all's TSCn-v^^nitrn ]) PDF with 
that of the single point iSpringel fc Farran ([2007) esti- 
mate. As noted in i i3.2.1l th ere is a large are a of param- 
eter space that the .Springel fc Farrad (|2007l ) result fails 
to represent. 

Similar to my analysis of the Bullet Cluster I perform 
the analysis with 2,000,000 Monte Carlo realizations. Pa- 
rameter estimates converge to better than a fraction of 
a percent with only 20,000 realizations. 

4.1. Musket Ball Observed System Properties 

I show the observed Musket Ball Cluster parameter 
PDF's in Figures |4l-|6l each th e result of analyses pre- 
sented by iDawson et al.l (|2012( ). I refer to their "south" 
subciuster as halo 1 and "north" subciuster as halo 2. 
The mass PDF's, Figure H are the result of an MCMC 
analysis where NFW halos were simultaneously fit to 
the weak lensing signal. The relative velocity distribu- 
tions ,_Jj£UTej5Lare_tlw result of a bootstrap error anal- 
ysis ([Beers et al.lll990t ) of the 38 and 35 spectroscopic 
members of the north and south subclusters, respectively. 
The projected subciuster separation distribution. Figure 
[HI is the result of a bootstrap error analysis of the re- 
cursively estimate d subclusters' galax y number density 
centroids (see e.g. iRandall et al.ll2008l for a description 
of this method). For each Monte Carlo realization in- 
dividual values are drawn randomly from each of these 
distributions. 

4.2. Musket Ball System Dynamics Results 

This more complete analysis confirms that the Mus- 
ket Ball Cluster merger is both significantly slower and 
further progressed compared to the Bullet Cluster, see 
Figure [T] To estimate a lower limit on how much further 
progressed I perform an additional Monte Carlo analysis 
for TSCq^^^^.^^ — rS'Cog^ii^t assuming the marginalized 
TSCo distributions (see Appendices iBl and [C[). This is 
a lower limit since the Musket Ball observations, unlike 
the Bullet Clusters observations, cannot rule out the case 
that its subclusters have reached the apoapsis and are on 
a return trajectory (61% of the realizations have TSCi 
less than the age of the Universe at z = 0.53). I find that 
the Musket Ball is at least O.Slg ^Gyr (3. 5^15 times) 
further progressed than the Bullet Cluster, see Figure 
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Fig. 5. — Relative radial subciuster velocity PDF's inferred 
from spectroscopic redshifts the Musket Ball Cluster galaxies 
HDawson et al.ll20T^ '). 
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Fig. 6. — Projected sep aration PDF of the M usket Ball subciuster 
galaxy density centroids HDawson et al.ir2012l ). 

[51 This is in line with th e more approximate 3-5 times 
estimate of IDawson et al.l ([2012D . The Musket BaU's rel- 
atively large TSCq means that it has potential for provid- 
ing tighter constraints on ctdm, since the expected offset 
between the galaxies and dark matter will increase with 
increasing TSCq. However the expected offset also in- 
creases as a function of the cluster mass surface density 
and collision velocity. Without running SIDM simula- 
tions it is difficult to know which of these parameters are 
most important. The complete Musket Ball Cluster pa- 
rameter estimates are summarized in Table[4|and plotted 
in Appendix [Cl 

Note that just as a temporal prior was justified for 
the Bullet Cluster based on the observed shock front 
and increased temperature/mass estimate, I could ap- 
ply a similar yet opposite prior to the Musket Ball since 
the temperature/mass estimate is consistent with the 
weak lensing inferred mass (additionally no shoc k front 
is observed) . According to [Randall et al.l ([2002D if the 
cluster mass and inferred X-ray temperature or lumi- 
nosity cluster mass are approximately the same then 
TSCq > 2tsc, which in the case of the Musket Ball means 
TSCq > 1.75 Gyr. While this is consistent with my 
TSCq estimate for the Musket Ball, it is not entirely 
appropriate to apply this prior since the X-ray observa- 
tions are relatively shallow and cannot c onfidently rule 
out a temperature and luminosity boost ([Dawson et al.l 
l2012f) . However it is conceivable that this line of reason- 
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Fig. 7. — The posterior of the Musket BaU Cluster's TSCo and 
''3D (*col ) parameters is shown in grayscale with dark and light blue 
contours representing 68% and 95% confidence, respectively. For 
comparison the gray dashed contours are the Bullet Cluster's 68% 
and 95% confidence intervals copied from Figure \3\ The Musket 
Ball Cluster occupies a much different region of merger phase than 
the Bullet Cluster, having both a slower relative collision velocity 
and being observed in a much later stage of merger. 




Fig. 8. — The histogram presents the TSCo^^^^^^t ~ '^■S'CoB„u,t 
distribution from random draws of the respective marginalized 
TSCo distributions; showing that the Musket Ball Cluster merger 
is at least 0.8^q 4Gyr (3.5^j 5 times) further progressed than the 
Bullet Cluster merger. The blac k dashed line is the biweight- 
statistic location I IBeers et al.1198"^ , the dark and light blue regions 
denote the bias-corrected 68% and 95% lower and upper confidence 
limits, respectively. 

ing would be applicable with deeper X-ray observations, 
either for the Musket Ball or similar dissociative mergers. 

5. SUMMARY AND DISCUSSION 

I have introduced a new method for determining the 
dynamic properties and associated uncertainty of dis- 
sociative cluster mergers given only the most general 
merger observables: mass of each subcluster, redshift of 
each subcluster, and projected separation the subclus- 
ters. I find that this method addresses the primary weak- 
nesses of existing methods, namely enabling accurate pa- 
rameter estimation and propagation of uncertainty near 
the collision state with a convergent solution achieved 
in ^6 CPU hours. I have confirmed that the two NFW 
halo model is capable of achieving the required percent 
level accuracy by direct comparison with N-body hydro- 
dynamic simulations. 

In applying this method to the Bullet Cluster I not only 



determined its merger dynamic parameters but found 
that the bulk of uncertainty in these parameters is due 
to uncertainty in a, the angle of the merger with respect 
to the plane of the sky. Analyses that fail to account for 
the uncertainty in a (all existing N-body simulations of 
the Bullet Cluster) will significantly underestimate the 
uncertainty in their results. This highlights the need to 
carefully select and model many possible realizations of 
the merger when trying to infer results from N-body sim- 
ulations of a real merger (I discuss this further in §5.1|) . 
I have also shown how ex post facto priors can easily 
be applied to the results of the default priors to further 
constrain the inferred dynamic properties. In particu- 
lar accurate measurement of the cluster gas properties 
can enable approximately a factor of two better con- 
straint on the dynamic properties of the merger, prin- 
cipally through added constraint on the time scale of the 
merger. 

I have also applied this method to the Musket 
Ball Cluster validating the approximate results of 
iDawson et al.l ()2012D . Comparing the dynamic proper- 
ties of the Musket Ball with those of the Bullet I have 
shown that the Musket Ball represents a significantly dif- 
ferent volume of merger phase space. The Musket Ball 
Cluster, being S.St^J further progressed than the Bul- 
let Cluster, could potentially provide tighter constraints 
on ctdm since the offset between galaxies and dark mat- 
ter should increase with time post-merger for (Tdm > 0. 
And the larger the expected offset, the better the dark 
matter constraint w hen applying a method similar to 
iRandall et all (f2008l) . 

5.1. Suggested Uses of Method 

While a general method for determining the dynamics 
properties of merging clusters has numerous applications, 
several are worth noting. As noted N-body simulations of 
specific merging clusters are computationally expensive; 
in particular one SIDM simulation of a single dissociative 
merger requires '-^I-IO million CPU hours (private com- 
munication, James Bullock). Thus it is currently unfea- 
sible to simulate all confirmed dissociative mergers. This 
method can be used to quickly determine which mergers 
provide the best ctdm constraining power, enabling an 
efficient use of limited computational resources. 

Additionally it is inappropriate to simply simulate one 
realization of a dissociative merger due to the broad 
range of merger phase space allowed by uncertainty in ob- 
served parameters, as discussed in detail in §3.2.11 Thus 
multiple simulations of each merger are required to prop- 
erly represent the allowed phase space. One could con- 
ceivably reduce the number of required simulations by 
using the results of this method to select representative 
merger realizations that uniformly sample the merger 
phase space of interest (e.g. cluster mass, 1130(^001), and 
TSCo); then weight the results of each simulated real- 
ization by the integral of the corresponding local phase 
space PDF, as determined by this method. For example, 
one could estimate the uncertainty distribution of the 
(Tdm constraint inferred from SIDM simulations of the 
Bullet Cluster by weighting the constraint from each re- 
alization, where a realization with ?;3D(icoi)=3100kms~^ 
and T5'Co=0.3 Gyr would receive greater weight than 
one with W3D(icoi)=4000kms-i and TS'Co=0.2 Gyr, see 
Figure |3l Thus the results of this method will not only 
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TABLE 4 

Musket Ball Cluster parameter estimates 



Parameter 


Units 


Location^ 


68% LCL-UCL*^ 


95% LCL-UCL*^ 


M2001 


-1/^14- TV T 

10 Mq 


3.2 


2.3 - 4.3 


1.6 - 5.5 


M20O2 


1 nl4 i\ /T 
lU Mq 


1.7 


1.1 — 2.4 


0.6 - 3.3 


Zl 




0.5339 


0.5333 - 0.5345 


0.5326 - 0.5352 


Z2 




0.5316 


0.5305 - 0.5324 


0.5294 - 0.5331 


'-^proj 


Mpc 


1.0 


0.9 - 1.1 


0.7 - 1.3 


a 


degree 


48 


28 - 67 


13 - 78 




Mpc 


1.6 


1.2 - 2.9 


0.9 - 5.5 




Mpc 


2.1 


1.5 - 3.8 


1.1 - 7.3 


f3D(tobs) 


km 


670 


390 - 1100 


140 - 1500 


f3D(tcol) 


km 


2300 


2000 - 2500 


1800 - 2800 


TSCo 


Gyr 


1.1 


0.7 - 2.4 


0.5 - 5.8 


TSCl" 


Gyr 


3.5 


2.0 - 7.2 


1.4 - 12.0 


T 


Gyr 


4.8 


2.9 - 10.4 


2.2 - 22.7 



^ Biweight-statistic location (see e.g. IBeers et al.lll990l ). 

^ Bias-correct ed lower and up per confidence limits, LCL and UCL respec- 
tively (see e.g. IBeers et al.lll99 01. 

61% of the realizations with a valid TSCq (i.e. less than the age of the Uni- 
verse at the cluster redshift) have a valid TSC'i , meaning that it is possible 
that the Musket Ball Cluster is being observed in the "incoming" state. 



inform efficient selection of realizations to model but will 
reduce the number of simulations required to properly 
sample the posterior PDF's. Nevertheless SIDM simula- 
tions of mock clusters need to be performed to determine 
how much acceptable values of ctdm affect the inferred 
merger dynamics properties. 

General merger dynamic properties are also impor- 
tant for understanding how cluster mergers relate to 
other physical phenomena, such as galaxy evolution and 
radio relics. It is well established that galaxy clus- 
ters play an important role in the evolution of their 
member galaxies, but it is still unclea r whether clus 
ter mergers trigger star formation (e.g. iMiller fc Owen 



2Uo aiOwen et all2005tlRirrari et al.lbooa lH wane fc Led 



20091 ) ■ quench it (iPoggianti et al.ll2Q04D . or have no im- 
mediate effect (|Qiung et al.l |2010( )~ Studying mergers 
at different TSC may resolve these seemingly conflict- 
ing results by discriminating between slow-working pro- 
cesses (e.g. galaxy harassment or strangulation) and fast- 
acting process (e.g. ram pressure stripping). Similarly, 
studying global merger dynamic properties may resolve 
the mystery of why many mergers have associated radio 
relics (e.g. Barrena e t al. 20 09; va n Weeren et al..,20li[) 
yet others don't (e.g. iRussell et al.ll2011i r 

5.2. Extensions to the Method 

While this method has advantages over existing meth- 
ods there is room for considerable improvement. For ex- 
ample the method could be improved through the elim- 
ination of some of the simplifying assumptions of the 
model (see §2.1|) . One could attempt to incorporate sub- 
cluster mass accretion phys i cs in a manor similar to the 
work of lAngus fc McGaughl ()2008[ ) or attempt to account 
for the possibility of a non-zero impact parameter. The 
latter is potentially more difficult to properly incorpo- 
rate. In addition to complicating the equations of motion 
it is unclear what impact parameter prior is reasonable. 
For example, given a two-subcluster merger it might be 
expected that near head-on collisions are more likely, but 
without a study of mergers in cosmological N-body sim- 
ulations it is unclear how much more likely. Additionally 



as the impact parameter increases the amount of gas dis- 
sociated from each subcluster will decrease, but without 
a systematic study of various mergers in hydrodynamic 
simulations it is unclear exactly what prior an observed 
large dissociation of gas should infer. 

Another significant extension to the model could be the 
inclusion of SIDM physics. As mentioned in the previous 
section, one of the promising uses of this method is to 
suggest which mergers might provide the best ctdm con- 
straining power. However one could take this a step fur- 
ther by including an analytic treatment of SIDM physics 
(e.g. iMarkevitch et all 120041 ) . thereby enabling analytic 
estimates of ctdm relevant effects for a given merger. 
Then this method could be used in conjunction with ob- 
served dissociative mergers to place direct constraints on 
(TDM- Due to the increased complexity of the physics in- 
volved it would be necessary to verify this extension with 
SIDM N-body simulations. 

Note: W. Dawson will make available, upon request. 
Python code implementing the discussed Monte Carlo 
method. 
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APPENDIX 

POTENTIAL ENERGY OF TWO TRUNCATED NEW HALOS 
Generically the potential energy of a two-halo system with center to center separation r is 

V{r)= /$i(r')dm2, 



(Al) 



where ^i{r') is the gravitational potential of halo 1 as a function of radial distance r' from th e cen ter of the halo 1 
to the m ass element of halo 2, dm2- I derive $i(r) for the case of a truncated NFW halo in i jA.ll The integral of 
equation I All can be approximated as a summation over N x N mass elements, m2-j, each with area dr x dO, where i 
and j range from TV — 1, 

N-lN-l 

where r^^ is the distance from the center of halo 1 to the 2""^ halo's mass element 77x2^^ , as derived in ijA.21 and e is 
the softening length which reduces the effects of artificial singularities. 

Truncated NFW Gravitational Potential 
For an axially symmetric mass distribution the potential can be expressed as a series of Legendre Polynomials 



2nG 



{n + l/2)r"+i 7o 



27rGr" 
n + 1/2 



r'^-"Pn(r')dr' 



where 



p„(r) = (n + 1/2) / p(r,6')P„ (cos 6*) sin 61 d6'. 



(A2) 
(A3) 



Assuming a spherical NFW halo 



PNFW[r) 



r/rs{l + rjrsY 



only the zero"^ order term of Equation I A3I remains 



and Equation I A2I reduces to 
$NFw(?') : 

*NFW(»^) ^ 



47rG r 



r 

A-nGp. 



P^Fwi.'r) = po{r) 

?''pNFw(''') dr' 

.12 

■ dr' + r 



r'/rsil + r'/rs)'- 



r'/rsil + r'/rsY 



■dr' 



Since I truncate the NFW halo at r2oo the 00 in the second integral becomes r2oo and 



NFWt 



(r) = 



GM2OO 



ln(l + r/r,)-^ 



+ I'200 



if < f20o; 
if r > r2oo- 



(A4) 



Mass Elements of a Truncated NFW Halo 
Given the differential mass elements for a spherically symmetric halo 

dm = 2np{r, 9)r'^ sm{9)d9 dr, 

and discretizing the mass into elements with lengths Sr — ^2002/-^ a-nd 60 = t:/N the halo 2 mass elements are given 

by 

/ p{r')r''^siii{9')d9' dr'. 



rriij = 27T 



i Sr 



For an NFW halo this becomes 

m^j = ^irpsrl [cos(j 59) - cos ((j + 1)59)] 



{i + l)5r 



i 5r 



In 



{i + l)5r + rs 



i 5r + 
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BULLET CLUSTER RESULT PLOTS 

This section contains the parameter results array plots for the Bullet Cluster case including the added temporal prior 
of tj3.2.2l For ease of display the parameters are grouped in three categories (Input, Geometry, and Velocity & Time) 
resulting in a six subplot results array, see Figure [9l The Input parameters consist of: M2001, ^2002) zi, Z2, and dproj, 
where halo 1 refers to the "main" subcluster and halo 2 refers to the "bullet" subcluster. The Geometry parameters 
consist of the randomly drawn a, and calculated c?3d, and c?max- The calculated Velocity & Time parameters consist 
of: TSCo, TSCi, and T. 



O 

o 



H 

o 
o 



Input 



Geometry Velocity & Time 



For ease of display the results array is divided into six subplots, Figures I10I414I The Input parameters consist of: M20O1 , 



Fig. 
M2002 

consist of: i>3D(iobs)> '!'3d(*coi)> TSCq, TSCi, and T. 



9. 

21, Z2, and rfproj- The calculated Geometry parameters consist of: 



and dmax- The calculated Velocity & Time parameters 
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Fig. 10. — Bullet Cluster marginalized Input vs. Input parameters result plots, for the case including the added temporal prior of il3.2.2l 
Dark an d light blue colors correspond to 68% and 95% confidence intervals, respectively. The black dashed line is the biweight-statistic 
location IIBeers et al.lll982l) . 
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Fig. 11. — Bullet Cluster marginalized Input vs. Geometry parameters result plots, for the case including the added temporal prior of 
i|3.2.2l Dark and light blue colors correspond to 68% and 95% confidence intervals, respectively. 
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F ig. 12 . — Bullet Cluster marginalized Input vs. Velocity & Time parameters result plots, for the case including the added temporal 
of ^3.2.21 Dark and light blue colors correspond to 68% and 95% confidence intervals, respectively. 
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a (degrees) d^D^Kbs) (Mpc) d^^^ (Mpc) 

F ig. 13 . — Bullet Cluster marginalized Geometry vs. Geometry parameters result plots, for the case including the added temporal prior 
of i|3.2.2l Dark a nd light blue color s correspond to 68% and 95% confidence intervals, respectively. The black dashed line is the biweight- 
statistic location IjBeers et al.lll98^ ). 
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Fig. 1 4. — B ullet Cluster marginalized Geometry vs. Velocity & Time parameters result plots, for the case including the added temporal 
prior of i|3.2.2l Dark and light blue colors correspond to 68% and 95% confidence intervals, respectively. 
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Fig. 15. — Bulle t Clus ter marginalized Velocity & Time vs. Velocity & Time parameters result plots, for the case including the added 
temporal prior of il3.2.2l Dark a nd light blue color s correspond to 68% and 95% confidence intervals, respectively. The black dashed line 
is the biweight-statistic location lIBeers et al.lll98^ . 
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MUSKET BALL CLUSTER RESULT PLOTS 

This section contains the parameter results array plots for the Musket Ball Cluster. Similar to f|B]the parameters are 
grouped in three categories {Input, Geometry, and Velocity & Time) resulting in a six subplot results array, see Figure 
[9l The Input parameters consist of: M2001, M2002; ^i; ^2, and dproj, where halo 1 refers to the "south" subcluster 
and halo 2 refers to the "north" subcluster. The calculated Geometry parameters consist of: a, dso, and dmax- The 
calculated Velocity & Time parameters consist of: f3D(iobs), ^^^{tcoi), TSCq, TSCi, and T. 




Fig. 16. — Musket Ball Cluster marginalized Input vs. Input parameters result plots. Dark an d light blue colo rs correspond to 68% and 
95% confidence intervals, respectively. The black dashed line is the biweight-statistic location (Beers et al.lll983) . 
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Fig. 17. — Musket Ball Cluster marginalized Input vs. Geometry parameters result plots. Dark and light blue colors correspond to 68% 
and 95% confidence intervals, respectively. 
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Fig. 18. — Musket Ball Cluster marginalized Input vs. Velocity & Time parameters result plots. Dark and light blue colors correspond 
to 68% and 95% confidence intervals, respectively. 
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Fig. 19. — Musket Ball Cluster marginalized Geometry vs. Geometry parameters result plots. Dark a nd light blue colo rs correspond to 
68% and 95% confidence intervals, respectively. The black dashed line is the biweight-statistic location ijBeers et al.lll98^ . 
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Fig. 20. — Musket Ball Cluster marginalized Geometry vs. Velocity & Time parameters result plots. Dark and light blue colors correspond 
to 68% and 95% confidence intervals, respectively. 
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Fig. 21. — Musket Ball Cluster marginalized Velocity & Time vs. Velocity & Time parameters result plots. Dark and light blue colo rs 
correspond to 68% and 95% confidence intervals, respectively. The black dashed line is the biweight-statistic location IjBeers et al.lll98^ . 
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